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ORIGINAL ARTICLE 

An Integrated Population Pharmacokinetic Meta-Analysis 
of Propofol in Morbidly Obese and Nonobese Adults, 
Adolescents, and Children 

J Diepstraten 1 , V Chidambaran 23 , S Sadhasivam 23 , HJ Blusse van Oud-Alblas 4 , T Inge 5 , B van Ramshorst 6 , EPA van Dongen 7 , 
AA Vinks 38 and CAJ Knibbe 1 ' 9 

This study describes a population pharmacokinetic meta-analysis of propofol to characterize the influence of body size 
measures and age in morbidly obese and nonobese adults, adolescents, and children. Sixty morbidly obese and nonobese 
adult patients (55-167 kg; 21-79 years) and 34 morbidly obese and nonobese adolescents and children (37-184 kg; 9-20 years) 
were included. The results show that clearance increased with total body weight in an allometric function while age was found 
to influence clearance in a bilinear fashion with two distinct slopes, reflecting an initial increase and subsequent decrease 
as a result of aging. Using these two functions, the influence of both (over)weight and age on propofol clearance was well 
characterized, which may provide a basis for dosing across this diverse group of patients. 
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BACKGROUND 

Although total body weight (TBW) of children and adolescents 
increases due to growth-related processes across childhood, 
obesity may also substantially contribute to increases in body 
weight. 1 As a result, morbidly obese children and adolescents 
may be as heavy as adults, even though growth-related pro- 
cesses have not yet been completed. The question then arises 
whether TBW, which is commonly used to adjust dosing in 
children and adolescents, is the appropriate measure to adjust 
doses of drugs in obese children and adolescents. Similarly for 
adults, there is a lively discussion about the best size descrip- 
tor for changes in pharmacokinetics due to obesity. 2 3 As little 
is known on how key pharmacokinetic parameters such as 
clearance change in morbidly obese children, adolescents, 
and adults as compared with their nonobese controls, studies 
are needed analyzing a wide range of ages and related TBWs. 

Propofol is widely used for induction and maintenance of 
anesthesia in both nonobese and (morbidly) obese adults, 
adolescents, and children. Recently, the pharmacokinet- 
ics of propofol have been compared in premature neonates 
and adults, 4 in morbidly obese and obese adults, 56 and in 
(morbidly) obese children and adolescents. 7 In all these stud- 
ies, TBW proved the most predictive covariate for clearance, 
either by using a standard allometric function 5-7 or a TBW- 
dependent exponent allometric function. 4 However, a meta- 
analysis on the basis of all data sets in morbidly obese adults, 
adolescents, and children together with their nonobese 
controls in which the influence of obesity and ageing is disen- 
tangled has not been performed. 



Therefore, the aim of this study was to perform a population 
pharmacokinetic meta-analysis of propofol combining data 
from morbidly obese and nonobese adults, adolescents, and 
children. To study how obesity and age influence pharma- 
cokinetic parameter estimates in this diverse patient group, 
specific emphasis was placed on the evaluation of the influ- 
ence of TBW, body mass index, ideal body weight, 8 lean body 
weight (LBW), 910 and/or age on the different pharmacokinetic 
parameters. 

RESULTS 
Subjects 

Ninety-four adults, adolescents, and children with a mean 
TBW of 94 kg (range: 37-1 84 kg) were included from which 
1,652 concentration measurements were available. Demo- 
graphic characteristics of the morbidly obese and nonobese 
patients are summarized in Table 1 . 

Pharmacokinetics 

A three-compartment pharmacokinetic model adequately 
described the time course of the propofol whole-blood concen- 
trations in all morbidly obese and nonobese adults, adoles- 
cents, and children. Exploratory plots of the tested covariates 
TBW, body mass index, ideal body weight, LBW, and age 
against individual post hoc parameter estimates of the sim- 
ple model without covariates (model A) showed a potential 
relation between clearance and TBW, with lower values for 
children and adolescents across the entire body weight range 
(Figure 1, model A). In addition, potential relationships were 
observed between central volume of distribution (V1) and 
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TBW or LBW, and between intercompartmental clearance 
from the central to the second peripheral compartment (Q3) 
and TBW (figures not shown). There were no visual trends 
between the explored covariates and other pharmacokinetic 
parameters in the simple model without covariates (model A). 

Subsequently, as depicted in Table 2, all body size mea- 
sures and age were separately incorporated on clearance, 
V1 , and Q3 in the model and tested for significance (see sec- 
tion Covariate analysis). The analysis showed that TBW was 
the most predictive covariate for propofol clearance when 
implemented using an allometric function (model B, decrease 
in objective function value (OFV) of 84.4 points; P < 0.001; 
Table 2). Figure 1 model B (and model A) shows that ado- 
lescents with the same TBW as adults had lower clearance 



values (gray vs. black symbols, respectively). Therefore, 
in model C, a separate value for propofol clearance in ado- 
lescents vs. adults was estimated. This resulted in another 
reduction in OFV by 23.5 points (P < 0.001) with individual 
clearance values for an adolescent of 70 kg and an adult of 
70 kg of 1.75ml/min and 2.18ml/min, respectively (Table 2, 
Model C). The nonlinear increase of propofol clearance 
with TBW proved the same for both groups and was best 
described using an allometric function with an estimated 
exponent of 0.73 (coefficients of variation percentage: 6.9) 
(Figure 1, Model C). 

However, when the simple model without covariates was 
evaluated for the effect of age (Figure 2, model A), it was 
found that clearance increased until the median age of 41 



Table 1 Baseline characteristics of all morbidly obese and nonobese adults, adolescents, and children included in the current analysis 



Adults Adolescents and children 





All Patients 


Morbidly obese 6 


Nonobese 11,12 


Obese 7 


Nonobese 13 


Number 


94 


20 


40 


20 


14 


Gender (M/F) 


30/64 


4/16 


16/24 


8/12 


2/12 


Age (years) 


38 (20) 


45(12) 


55(12) 


16(2) 


14(3) 


TBW (kg) 


94 (35) 


124 (20) 


74(11) 


125 (29) 


54(13) 


BMI (kg/m 2 ) 


33(12) 


43 (6) 


26 (4) 


46 (9) 


21 (6) 


IBW (kg) 


61 (9) 


61 (7) 


64 (8) 


59(12) 


55 (9) 


LBW (kg) 


54(14) 


60 (9) 


50(10) 


63(14) 


37 (8) 



Data are presented as mean with standard deviation (SD). 

BMI, body mass index; F, female; IBW, ideal body weight; 8 LBW, lean body weight; 9 M, male; SD, standard deviation; TBW, total body weight. 



Model A Model B 




Total body weight (kg) Total body weight (kg) 



Figure 1 Individual post hoc propofol clearance estimates vs. total body weight for the simple model (model A) and three covariate 
pharmacokinetic models (B, C, and D) for morbidly obese adults (black circles), adolescents and children (gray circles) and their nonobese 
controls (n = 94). In model B, the black line indicates the population clearance values for both the adult and adolescent population; in model 
C, the black line indicates the population clearance values for adults and the gray line the population clearance values for adolescents; and in 
model D, the black dotted lines indicate the population clearance values for 15, 41 , and 65 years. 
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Table 2 Results of covariate analysis for the three-compartment pharmacokinetic model of propofol in the combined data set of morbidly obese and nonobese 
adults, adolescents, and children 



Parameter 


Model 


Model description 


Relationship of covariate 


Number of structural parameters 


A OCX/ 

AUrv 




A 


Simple model 




11 




UL 




Age bilinear 


CL. = CL x F a 

i pop age 


13 


-16.7 


UL 




Lbvv (yj linear 


w /i n\A/ /cf>\ 
UL. = UL pop x (LBW/60) 


1 1 


-45.8 


CL 




1 bw linear 


CL. = CL pop x (TBW./70) 


1 1 


-68.5 


CL 


B 


TBW allometnc 


.... /tdiii nr\\7. 

CL. = CL pop x (TBW/70) Z 


12 


-84.4 


CL 


C 


TBW allometric, CL adolescents , and CL adults 


CL. = CL pop ado|escents x (TBW/70) Z 

CL i = CL pop, adults x(TBW/70) z 


13 


-107.9 


CL 


D 


TBW allometric and age bilinear 


CL = CL pop x(TBW/70) z xF age a 


14 


-116.1 


V1 




TBW linear 


V1 i = V1 pop x(TBW/70) 


11 


-0.2 


V1 




LBW 9 linear 


V1 i = V1 pop x (LBW/60) 


11 


-5.8 


Q3 




TBW linear 


Q3 i = Q3 pop x(TBW/70) 


11 


-18.1 


Q3 




TBW allometric 


Qa = Q3 pop x(TBW/70) d 


12 


-18.1 


Final model 


E 


TBW allometric and age bilinear on CL 


CL = CL pop x(TBW/70) z xF age * 


14 


-143.0 


CL and Q3 




and TBW linear on Q3 








Q3 i = Q3 pop x(TBW/70) 







a Age < 41 years: F age = (1 + b x (Age - 41)) and Age > 41 years: F age = (1 + c x (Age - 41)). 

CL, clearance; CL, clearance in /th individual; CL pop , population mean value for clearance; F age , age factor for clearance; LBW, lean body weight; AOFV, delta objec- 
tive function value as compared with simple model; Q3, compartmental clearance between V1 and V3; TBW, total body weight; V1 , central volume of distribution; 
z, allometric scaling factor in CL = CL 70kg x (TBW/70) Z . 
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Figure 2 Interindividual variability of propofol clearance vs. age for 
the simple model without covariates (model A) and the final covariate 
model including age and total body weight on propofol clearance 
(model E). 



years after which it decreased. As a result, instead of esti- 
mation of two different population values for adolescents vs. 
adults as in model C, in model D, age was implemented using 
a bilinear function which significantly reduced the OFV (AOFV 
as compared with model C = -8.2 points; P < 0.005). On the 
basis of the covariates of model D, the interindividual vari- 
ability of propofol clearance was reduced by 50%. Figure 2 
model E shows that after implementation of age in a bilinear 
function, interindividual variability was randomly distributed 
with age. Figure 1, model D shows the post hoc propofol 
clearance estimates for model D vs. TBW with population 
predictions for clearance for three different ages (15, 41 , and 
65 years), illustrating the bilinear relation with age in model D. 
The final equation for propofol clearance was given as follows 
(Eq.1): 

CL i= CL pop x (TBW/70)° 77 x F age a (1) 
a Age < 41 years: F age = (1+ 0.0103 x (Age - 41)) 
Age > 41 years : F age = (1 - 0.00539 x (Age - 41)) 

where CL represents clearance in the /th individual, CL pop is 
the population mean value for clearance in an individual of 
70 kg and of 41 years, TBW i is the TBW of the /th individual, 
and 70 is the standard body weight in kilograms. 

Concerning covariates for V1 , Table 2 shows that there was 
only a modest influence of the body size descriptors on V1 
with a trend toward an increase in V1 with LBW (P > 0.005). 
There was a substantial shrinkage (43%) on V1, which not 
only renders plots using post hoc parameter estimates less 
reliable but also indicates that the individual data in the data 
sets are not rich in information about this parameter. 11 There- 
fore, no covariate on V1 was incorporated in the final model. 
By contrast, TBW as covariate for Q3 significantly improved 
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the model (AOFV = -1 8.1 ; P < 0.005; Table 2) and was there- 
fore considered the final model (Table 2, model E). There was 
no influence of the explored covariates on the other pharma- 
cokinetic parameters (Q2 and V2). 

Table 3 lists all parameter estimates including their coef- 
ficient of variation values and OFVs of the simple model 
(model A) and the final model (model E). The observed vs. 
population-predicted plots stratified by the different cohorts 
in Figure 3 confirm that the final model describes not only 
the study population as a whole but also the individual study 
populations without bias. The stability of the final model was 
shown by the bootstrap analysis (Table 3). 

Figure 4 shows population propofol clearance values vs. 
age for different TBWs using the final model (model E). This 
figure shows both the allometric increase of propofol clear- 
ance with TBW as the distance between the weight classes 
decreases with increasing TBW, and the bilinear relationship 
of propofol clearance with age. 

DISCUSSION 

To describe the influence of obesity and age on the phar- 
macokinetics of propofol, a population pharmacokinetic 



meta-analysis was performed using data from morbidly 
obese adults, adolescents, and children and their nonobese 
controls. In this study, a wide range of TBW (37-1 84 kg) and 
age (9-79 years) was studied, with data from (morbidly) 
obese and nonobese individuals in each age range. The 
results of the systematic analysis shows that a combination 
of TBW and age proved to best capture changes in propo- 
fol clearance as a result of obesity and ageing. Although it 
is yet unknown how these results should be put in physi- 
ological perspective, the current model seems to provide 
the best description of the data from these largely divergent 
patient populations. 

In recent reports in (morbidly) obese adults, it was 
shown that the increase in propofol clearance was related 
to TBW and could be best described using an allometric 
function. 56 In addition, an allometric relationship between 
TBW and propofol clearance was found in a data set of 
morbidly obese adolescents. 7 Allometric scaling factors of 
0.72 (ref. 6) and 0.80 (ref. 7) were estimated for morbidly 
obese adults and children and adolescents, respectively. 
As these factors are close to the factor of 0.75 predicted 
by allometry theory, 12 this implies that obese individuals 
can be viewed as "large individuals" (a different body size) 



Table 3 Population pharmacokinetic parameter estimates for the simple and the final pharmacokinetic model for propofol in nonobese and (morbidly) obese 
children, adolescents, and adults including the bootstrap values of the final model 




Simple model 


Final model 


Bootstrap final model 


Parameter 


Mean (CV%) 


Mean (CV%) 


Mean (CV%) 


Model 


A 


E 




Number of patients 


94 


94 


250 


CL (l/min) 


2.37 (3.8) 






CL 70 k g, 41 years ( l/min ) 3 




2.34 (4.3) 


2.31 (4.6) 


Z 




0.77 (6.9) 


0.77 (7.4) 


tiP 




0.0103(13.5) 


0.0094(18.5) 


& 




-0.00539 (-33.8) 


-0.00485 (-39.7) 


V1 (L) 


3.33(11.5) 


3.17(11.3) 


3.16(10.0) 


V2 (L) 


6.55(14.5) 


5.89(15.0) 


5.78(15.3) 


V3 (L) 


118(9.1) 


116(7.5) 


117(6.7) 


Q2 (l/min) 


1 .74 (9.7) 


1.60(11.7) 


1.57(12.2) 


Q3 (l/min) 


2.00(7.6) 






Q3 70kg (l/min)' 




1.50 (6.2) 


1.46 (5.0) 


Interindividual variability (%) 








CL 


35.1 (15.5) 


17.5(13.9) 


17.7(12.6) 


V1 


46.6 (44.4) 


50.6 (41.3) 


51.1 (38.8) 


V3 


40.4 (42.0) 


36.2 (34.9) 


35.4 (27.8) 


Q3 


48.1 (42.3) 


40.4 (37.5) 


35.1 (34.4) 


Proportional intraindividual error (%) 


24.3(10.4) 


24.3(10.3) 


24.1 (11.1) 


OFV 


-2,331 


-2,474 


-2,500 



aCL = CL 70kg x(TBW/70)*xF ao 



: (1 + b x (age - 41 )) and Age > 41 y: F = (1 + c x (age - 41 )). 



b Age < 41 y: F ac 
C Q3.- 

b, age factor for clearance age < 41 years; c, age factor for clearance age : 



of 70 kg; CL 



41 years; CL, clearance; CL 70kg , population mean value for clearance in an individual 



, population mean value for clearance in an individual of 70 kg and 41 years; CV, coefficient of variation of the parameter values; F a age 



factor for clearance; OFV, objective function value; Q2, compartmental clearance between V1 and V2; Q3, compartmental clearance between V1 and V3; Q3 70kg , 
population mean value for compartmental clearance between V1 and V3 in an individual of 70 kg; V1 , central volume of distribution; V2, peripheral volume 1 ; V3, 
peripheral volume 2; z, allometric scaling factor in CL. = CL x (TBW/70) Z . 
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Figure 3 Observed vs. population-predicted In propofol concentrations of the final model (model E). Panels represent data of morbidly obese 
adults, nonobese adults, morbidly obese children and adolescents, and nonobese children and adolescents. The solid gray line represents 
the line of identity, x = y. 

instead of individuals "having excess body fat" (a differ- 
ent body composition). 2 Although these results were con- 
firmed in the current meta-pharmacokinetic analysis, we 
also showed that morbidly obese adolescents cannot be 
viewed as "adults" as their propofol clearance proved lower 
than that of morbidly obese adults with the same TBW 
(Figure 1, model A). This difference in propofol clearance 
could be described with two separate functions for propofol 
clearance; i.e., one equation for children and adolescents 
and one equation for adults (model C). Alternatively and 
significantly better, age was incorporated as covariate on 
propofol clearance using a bilinear function (models D and 
E). Therefore, in the final model, the influence of age and 
obesity on propofol clearance was described using both 
TBW and age as covariates for propofol clearance. This 
final equation (Eq. 1) is independent of the definitions for 
age (e.g., adolescents and adults) and obesity (e.g., obese 
and morbidly obese) categories and might prove useful for 
clinical practice. 

In this study, there was no significant relationship 
between body size measures and volumes of distribution. 
Previously, age and TBW have been identified as covari- 
ates for volumes of distribution of propofol in nonobese 
and obese patients. 5 ' 1314 As a result of the finding that LBW 
correlated with V1, Ingrande etal? 5 suggested to use LBW 
for the induction of anesthesia with propofol. The lack of 
significant influence of LBW on volume of distribution in 
our analysis may be explained by the fact that the studies 
included in the current analysis mainly contained observa- 
tions following propofol maintenance infusions. As such, 



20 



40 

Age (years) 



60 



80 



Figure 4 Model-based predictions of population clearance estimates 
of propofol vs. age for patients with different total body weights. 

these data sets did not contain sufficient observations just 
after the induction bolus dose of propofol to adequately 
describe early (re-)distribution and the influence of covari- 
ates on volumes of distribution. It therefore seems that 
additional research is needed to characterize covariates 
predictive of volume of distribution that will allow estima- 
tion of propofol-loading doses in morbidly obese adults and 
children. 

This study had a few limitations. We investigated a cohort 
of children and adolescents with a lower age limit of 9 
years. In clinical practice, propofol is administrated to even 
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younger children, and therefore, more research is needed 
to investigate if the current findings are also applicable for 
children younger than 9 years. In addition, the pharmacody- 
namics of propofol were not included in the current analy- 
sis. Obesity influenced the pharmacodynamics of propofol 
in children and adolescents using propofol. 16 Underlying 
diseases such as diabetes and a changed (pathophysi- 
ology in (morbidly) obese patients may also influence the 
pharmacodynamics of propofol. 17 To develop dosing algo- 
rithm for propofol in morbidly obese adults, children, and 
adolescents, an integrated pharmokinetic and pharmaco- 
dynamic meta-analysis is urgently needed. 

It remains to be speculated how the influence of TBW 
on propofol clearance that was found in our study can be 
explained. Studies have shown that obese patients suffer 
from low-grade inflammation, 18 which is probably the under- 
lying cause of the high prevalence of nonalcoholic steato- 
hepatitis. 19 It is known that nonalcoholic steatohepatitis 
increases fat deposition in the liver causing sinusoidal nar- 
rowing and altered functional morphology of the liver. 20 By 
contrast, because of increased blood volume and cardiac 
output, hepatic blood flow is possibly increased in obese 
subjects. 21 As a result, increased propofol clearance may 
be anticipated as propofol is a high extraction ratio drug 22 
mainly metabolized by various UDP-glucuronosyltransfer- 
ase enzymes. 23 Data on other high-extraction drugs and 
drugs metabolized by UDP-glucuronosyltransferase sug- 
gest that both UDP-glucuronosyltransferase activity 24-26 
and liver blood flow 2728 are increased in obese adults. 
Furthermore, UDP-glucuronosyltransferase activity is 
increased in obese adolescents as compared with non- 
obese adolescents. 29 Even though this cannot be proven, 
it can be hypothesized that hepatic blood flow is even more 
increased due to prolonged duration of obesity in adults 
as compared with adolescents with the same TBW. This is 
supported by the fact that age could be incorporated as a 
covariate on propofol clearance. As propofol clearance is 
limited by the blood flow through the liver, the effect of both 
TBW and age on propofol may be explained by changes in 
liver blood flow. 

In this pharmacokinetic meta-analysis, we developed a 
model for scaling propofol clearance over wide ranges of 
TBW and age using data from morbidly obese adults, ado- 
lescents, and children and their nonobese controls. The 
results show that TBW was the most predictive covariate for 
propofol clearance among patients when implemented as an 
allometric function. In addition, age was incorporated using a 
bilinear function with two distinct slopes, reflecting an initial 
increase and subsequent decrease in clearance as a result 
of age. Using these two functions, the influence of both (over) 
weight and age on propofol clearance was well character- 
ized, which may provide a basis for dosing across this diverse 
group of patients. 

METHODS 

Patients. Data of five previously published studies were used 
for this analysis. 6 ' 730-32 Patient characteristics of the five 



different studies are provided in Table 1 . Details of the studies 
are briefly summarized when relevant to the current analysis. 

Morbidly obese adults. Twenty morbidly obese adults sched- 
uled for bariatric surgery with a mean TBW of 1 24 kg (range: 
98-1 67 kg) received either a propofol induction dose of 200 
or 350 mg. Maintenance propofol infusion rate was initiated 
at 10mg/kg times TBW/h and adjusted to keep Bispectral 
index values between 40 and 60. Remifentanil continuous 
infusion was administrated 25 ug/h/kg based on ideal body 
weight. Multiple blood samples were collected before the 
start of the propofol bolus until 150min after the end of the 
infusion. 6 

Nonobese adults. Forty nonobese adults with a mean TBW 
of 74 kg (range: 55-98 kg) were included. Twenty-four female 
patients received a bolus injection of 2.5mg/kg of propofol 
for induction of anesthesia while anesthesia was maintained 
with isoflurane. Following onset of unconsciousness, fentanyl 
(0.003-0.005 mg/kg) was administrated and as required addi- 
tional amounts of 50-100 mg. Blood samples were collected 
from 1 min until 180min after the induction dose of propofol. 30 
Of these 24 patients, only 20 patients were in this study 
because height measure of 4 patients was not available. 
Another 20 nonobese intensive care patients received con- 
tinuous propofol infusions for 2-5 days with propofol doses 
based on the Ramsay six-point-scale morphine administra- 
tion as an analgesic was given when the patient was con- 
sidered to be in pain (seven patients; each patient received 
morphine as a continuous infusion of 40-50 mg/day). Blood 
samples were collected four times daily during propofol main- 
tenance infusion for 2-5 days. 31 

Morbidly obese children and adolescents. In 20 morbidly 
obese adolescents and children scheduled for bariatric 
surgery with mean TBW of 125 kg (range: 70-1 84 kg) and 
mean age of 16 years (range: 9-18 years), propofol was 
administered using dosing weight calculated according to 
the method of Servin et al. Fentanyl 100 ug was adminis- 
tered just after the induction, and 50 ug doses were admin- 
istered in case of inadequate analgesia. Blood samples 
were collected before the start of the propofol infusion and 
from 1 5 min after the start of infusion until 1 20 min after the 
end of the infusion. 7 

Nonobese children and adolescents. In 14 nonobese ado- 
lescents and children with mean TBW of 54 kg (range: 
37-82 kg) and a mean age of 14 years (range: 9-20 
years), anesthesia was induced with a bolus dose of pro- 
pofol (4 mg/kg) and maintained with propofol by continu- 
ous infusion (2-1 0 mg/kg/h) and remifentanil (0.2-1 ug/ 
kg/min) for scoliosis surgery. Blood samples were taken 
before induction of anesthesia, and at -15 or 30 min after 
the start of propofol infusion until 120 min after the end of 
the infusion. 32 

Data analysis and internal validation. The analysis was per- 
formed by means of nonlinear mixed-effects modeling using 
NONMEM (version VI, release 1.1; GloboMax LLC, Hanover, 
MD) 33 with S-plus (version 6.2; Insightful Software, Seattle, 
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WA) to visualize the data. Discrimination between different 
models was made by comparison of the OFV (-2 log likeli- 
hood (OVF)). A value of P < 0.05, representing a decrease 
of 3.8 in the OVF, was considered statistically significant. 
In addition, goodness-of-fit plots (observed vs. individually 
predicted, observed vs. population predicted, conditional 
weighted residuals vs. time, and conditional weighted residu- 
als vs. population predictions) were used for diagnostic pur- 
poses. Furthermore, the confidence interval of the parameter 
estimates, the correlation matrix, and visual improvement 
of the individual plots were used to evaluate the model. 
r|-shrinkage as defined by Karlsson etal. 34 was calculated for 
all model parameters for which interindividual variability was 
estimated. The internal validity of the population pharmaco- 
kinetic models was assessed by a per-study stratified boot- 
strap resampling method using 250 replicates. 33 Validation 
using external data sets was not part of the current analysis. 

Pharmacokinetic model. Log-transformed propofol concen- 
tration data were described by a three-compartment model 
(NONMEM VI, ADVAN11, and TRANS4) parameterized in 
terms of volume of distribution of the central (V1), first (V2), 
and second peripheral compartments (V3), intercompart- 
mental clearance from the central to the first peripheral com- 
partment (Q2) and from the central to the second peripheral 
compartment (Q3), and clearance from the central compart- 
ment (CL). 

The interindividual value (post hoc value) of the parameters 
of the /th subject was modeled by the following formula: 



ft = 0 m 



(2) 



where 6 is the population mean and ri. is a random vari- 

mean 1 1 'i 

able with mean zero and variance co 2 , assuming lognormal 
distribution in the population. 

The intraindividual variability, resulting from assay errors, 
model misspecifications, and other unexplained sources, was 
best described with a proportional error model. This means 
for the /th observed log-transformed propofol concentration of 
the /th individual, the relation (Y): 



n='ogc pi 



(3) 



where c pred is the predicted propofol concentration and e is 
the random variable with mean zero and variance a 2 . 

Covariate analysis. Covariates were plotted independently 
against the individual post hoc parameter estimates of all 
pharmacokinetic parameters and the conditioned weighted 
residuals to visualize potential relations. The following covari- 
ates were tested: TBW, body mass index, ideal body weight 8 
and LBW, 9 ' 10 gender, and age. Covariates were tested using 
linear and power equations: 



P = P„ 



Cov 



Cov c1 



(4) 



where P. and P p represent individual and population 
parameter estimates, respectively, Cov represents the 



covariate, and Cov standard represents a standardized (i.e., 
70 kg for TBW) or median value of the covariate for the 
population, z represents the scaling factor, which was fixed 
to 1 for a linear function or an estimated value for a power 
equation. 

The influence of the covariate age on clearance was also 
tested using a bilinear function with two distinct slopes, i.e., a 
linear increase in clearance for age values below the median 
age and a linear decrease in clearance for age values higher 
than the median age 35 (Eq. 5). 



CL- CL pop x F age 
F age (age < median age) 

F age (age > median age) = 



(5) 



1+b x 



1+ C X 



age (in years) 
- median age (in years) 
age (in years) 
- median age (in years 



Potential covariates were separately entered into the model 
and statistically tested by use of the OFV and if applicable 
the 95% confidence interval of the additional parameter. A 
P < 0.005 was applied to evaluate the covariates in the for- 
ward inclusion (OFV decrease > 7.9), whereas the backward 
deletion procedure used a stricter criterion (OFV decrease > 
10.8; P < 0.001). When more than one significant covariate 
for the simple model was found, the covariate-adjusted model 
with the largest decrease in objection function was chosen 
as a basis to sequentially explore the influence of additional 
covariates using the same criteria. Finally, after forward inclu- 
sion, a backward exclusion procedure was applied to justify 
the covariate. The choice of the covariate model was further 
evaluated as discussed under the section Data analysis and 
internal validation. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

S There is very limited knowledge on how key 
pharmacokinetic parameters such as clear- 
ance change in morbidly obese children, 
adolescents, and adults as compared with their 
nonobese counterparts. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

S This study characterizes the influence of both 
morbid obesity and age on the pharmacokinet- 
ics of propofol using data from patients varying 
in TBW from 37 to 1 84 kg and in age from 9 to 
79 years. 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y Propofol clearance increased with TBW in an 
allometric function, whereas age was found to 
influence clearance in a bilinear fashion with 
two distinct slopes, reflecting an initial increase 
and subsequent decrease as a result of aging. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y While the proposed pharmacokinetic model 
was the best description of the available data, 
the results can be used to further explore 
physiological explanations for the findings. 
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